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Abstract 

We introduce a new technique for following the formation and evolution of galaxies 
in cosmological N-body simulations. Dissipationless simulations are used to track the 
formation and merging of dark matter halos as a function of redshift. Simple prescrip- 
tions, taken directly from semi-analytic models of galaxy formation, are adopted for 
gas cooling, star formation, supernova feedback and the merging of galaxies within the 
halos. This scheme enables us to explore the clustering properties of galaxies and to 
investigate how selection by luminosity, colour or type influences the results. In this 
paper, we study the properties of the galaxy distribution at z = 0. These include B and 
K-band luminosity functions, two-point correlation functions, pairwise peculiar veloci- 
ties, cluster mass-to-light ratios, B — V colours and star formation rates. We focus on 
two variants of a cold dark matter (CDM) cosmology: a high-density ($7 = 1) model 
with shape-parameter T = 0.21 (rCDM), and a low-density model with 0, = 0.3 and 
A = 0.7 (ACDM). Both models are normalized to reproduce the I-band Tully-Fisher 
relation of Giovanelli et al. (1997) near a circular velocity of 220 km s^^. Our results 
depend strongly both on this normalization and on the adopted prescriptions for star 
formation and feedback. Very different assumptions are required to obtain an acceptable 
model in the two cases. For rCDM, efficient feedback is required to suppress the growth 
of galaxies, particularly in low- mass field haloes. Without it, there are too many galax- 
ies and the correlation function exhibits a strong turnover on scales below 1 Mpc. For 
ACDM, feedback must be weaker, otherwise, too few L^, galaxies are produced and the 
correlation function is too steep. Although neither model is perfect both come close to 
reproducing most of the data. Given the uncertainties in modelling some of the critical 
physical processes, we conclude that it is not yet possible to draw firm conclusions about 
the values of cosmological parameters from studies of this kind. Further observational 
work on global star formation and feedback effects is required to narrow the range of 
possibilities. 
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1 Introduction 



A major motivation for carrying out N-body simulations of large-scale structure formation 
is to test theories for the origin of structure and to estimate cosmological parameters, such 
as the density fl, or the cosmological constant A. Physically accurate identification of the 
positions, velocities and intrinsic properties of galaxies is necessary if such simulations are to 
provide estimates of statistics such as the spatial and velocity correlations of galaxies, which 
can be reliably compared with observations. 

Dissipationless simulations of gravitational clustering using tens of millions of particles are 
now carried out routinely on parallel supercomputers. Such simulations are able to resolve 
the formation and evolution of the dark matter halos of typical galaxies over cosmologically 
significant volumes. However, once these galaxy-sized halos merge to form larger structures, 
such as groups or clusters, they are quickly disrupted and are no longer distinguishable as 
separate entities within the more massive systems. This is commonly referred to as the 
"overmerging problem". Several recent papers have demonstrated that with sufficient force 
and mass resolution (~ 1 — 3 kpc and 10*^ — 10^ Mq) the central cores of many simulated galaxy 
halos do "survive" in groups and clusters, but even with this resolution, substructure is still 
erased in the dense central regions (Tormen, Diaferio & Syer 1998; Ghigna et al 1998; Klypin 
et al. 1998). 

Dissipationless simulations do not address the fact that on scales of a few kiloparsecs, 
gas-dynamical processes are known to play a key role in determining the structure of galax- 
ies. Numerical simulations that treat the physics of the baryonic component can be used 
to overcome this deficiency. Indeed, it has been demonstrated that when radiative cooling is 
included, many individual dense knots of cold gas do survive within a present-day cluster (e.g. 
Carlberg, Couchman & Thomas 1990; Katz, Hernquist & Weinberg 1992; Navarro & White 
1994; Evrard, Summers & Davis 1994; Frenk et al. 1996). Such simulations require a heavy 
investment of CPU time and with present-day computer technology, it is difficult to follow 
the formation of galaxies in volumes that are large enough to study large-scale structure. In 
addition, the galaxies themselves are very poorly resolved and processes such as star formation 
and supernova feedback still have to be put into the simulations "by hand" because they are 
strongly influenced by structure much smaller than the resolution limit. 

In this paper, we introduce a technique for following the formation and evolution of galaxies 
in cosmological N-body simulations. Dissipationless simulations are used to track the forma- 
tion and merging of dark matter halos as a function of redshift. The most bound particle in 
each halo is identified as the site where cold gas condenses and forms stars. Later, when two 
or more halos merge in the simulation, these particles maintain their identities as separate 
galaxies except they can merge with the new central object on a timescale set by dynamical 
friction. Simple prescriptions are adopted for gas cooling rates, star formation and supernova 
feedback. These are based on simplified models of the physical processes or on simulation re- 
sults and are are taken directly from earlier semi-analytic studies of galaxy formation. These 
used analytic models rather than N-body simulations to specify the merging histories of dark 
matter halos for a given set of cosmological initial conditions. Such semi-analytic models were 
able to account for many aspects of the present-day galaxy population, for example the lumi- 
nosities, colours and morphologies of galaxies and the observed correlation of these properties 
with environment (White & Frenk 1991; Lacey & Silk 1991; Kauffmann, White & Guiderdoni 
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1993; Cole et al 1994; Baugh, Cole & Frcnk 1996b; Somerville & Primack 1998). In addi- 
tion, semi-analytic models have been used to study the evolution of the galaxy population to 
high redshift and to make predictions for faint galaxy counts, for the redshift distributions 
of faint galaxies, for the evolution of the morphology-density relation and for the abundance 
and properties of galaxies and damped Lya absorbers at high redshifts ( Lacey et al. 1993; 
Kauffmann, Guiderdoni & White 1994; Heyl et al 1995; Kauffmann 1995a,b; Baugh, Cole & 
Frenk 1996a; Kauffmann 1996a,b; Kauffmann & Chariot 1998; Baugh et al 1998; Mo, Mao & 
White 1998). 

The combination of semi-analytic galaxy formation models with cosmological N-body sim- 
ulations results in an important advance. It becomes possible to track the positions and 
velocities of galaxies as a function of time, in addition to properties such as stellar mass, 
luminosity, colour and morphology. One can then compute the standard statistical measures 
used to quantify the clustering of galaxies, such as the two-point correlation function or the 
pairwise peculiar velocities, and investigate how the selection of galaxies by luminosity, colour 
or type influences the results. One can extract mock catalogues from the simulations in order 
to study the clustering properties of galaxies in redshift space. One can investigate whether 
the galaxy distribution is biased relative to that of the dark matter and how this bias affects 
attempts to measure the density parameter Q from galaxy surveys. Finally, one can also study 
how the clustering of galaxies evolves with redshift. 

Some, but not all of these issues have been explored in previous papers. White et al (1987) 
assumed that a galaxy formed at the centre of a halo during its initial collapse and then tracked 
the positions of the central particles as a function of time. The circular velocities of galaxies 
were given by the the circular velocities of the halos in which they formed. Galaxies were 
also able to merge if they approached within a certain critical radius of each other. White et 
al (1987) demonstrated that at the present day, galaxies with high circular velocities cluster 
more strongly than galaxies with low circular velocities. Van Kampen & Katgert (1997) used 
a similar approach to follow the formation of galaxies within N-body simulations of cluster 
formation. In their scheme, each group of bound particles in virial equilibrium was replaced by 
a single massive particle with a softening corresponding to the radius of the group. The studies 
of White et al and Van Kampen & Katgert did not include any scheme for star formation and 
no attempt was made to compute the luminosities or colours of the galaxies in the simulation. 

Modelling of cooling, star formation, feedback and stellar evolution is included in the 
cosmological Eulerian grid calculations of Gen & Ostriker (1992,1993). Even though these 
calculations have insufficient resolution to follow the collapse and dynamics of galaxy halos, 
and so cannot reliably identify individual galaxies, Cen & Ostriker were able to address some 
of the issues we discuss in the present paper (e.g. bias and the dependence of clustering on 
galaxy age). A number of their conclusions prefigure our own. 

The results of semi-analytic models have been used to place galaxies within individual 
outputs of N-body simulations with higher resolution. Analytic methods are used to compute 
the luminosity function of galaxies within dark matter halos of given mass M and galaxies are 
then assigned randomly to particles in each simulated halo of this mass. Kauffmann, Nusser 
& Steinmetz (1997) used this method to study present-day bias as a function of galaxy lumi- 
nosity, colour and type, while Governato et al. (1998) used the same technique to analyze the 
clustering properties of Lyman break galaxies at redshifts ~ 3. The disadvantage is that the 
properties of the galaxies within each simulated halo do not depend on the merging history of 
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that particular halo. Correct results are only obtained if an average is taken over many halos. 
Furthermore there is no detailed correspondence between the galaxy populations assigned to 
the simulations at different times so that many questions about the evolution of galaxy clus- 
tering cannot be directly addressed. Finally, the method is subject to all the inaccuracies of 
the Press-Schechter theory and its extensions (see section 5 below). It is therefore advanta- 
geous to compute the merging histories of dark halos directly from the simulations. This has 
been done by Roukema at al. (1997) who studied the evolution of both halos and galaxies 
in simulations with scale-free initial conditions. Their method included prescriptions for star 
formation and they were able to demonstrate that different assumptions could strongly affect 
the predicted number density of low-luminosity galaxies. 

This paper is the first in a series studying the properties of galaxies identified in N-body 
simulations with cold dark matter (CDM) initial conditions. These simulations are much 
larger than those analyzed in previous work and are able to resolve the detailed merging 
history of the dark matter halo of an L^, galaxy within a volume of ~ 10^ Mpc^. We first 
describe the simulations and the methods used to construct halo merging trees. We outline 
the prescriptions adopted to model cooling, star formation, supernova feedback, galaxy-galaxy 
merging and the evolution of the stellar populations within galaxies. We show how the lumi- 
nosity functions of galaxies in the simulations compare with those derived from the analytic 
merging trees of Kauffmann & White (1993). 

We then study the properties of the galaxy distribution at z — 0. These include the I-band 
Tully-Fisher relation, B and K-band luminosity functions, the two-point spatial correlation 
function ^(r), two-point velocity correlations, cluster M/L ratios, the B — V colour distribu- 
tions of galaxies of different types, and the predicted Ha luminosity function. We focus on 
two variants of a cold dark matter (CDM) cosmology: a high-density {Q = 1) model with 
shape-parameter F = 0.21 (rCDM), and a low-density model with Q — 0.3 and A = 0.7 
(ACDM). The normalization (jg is chosen to match the observed abundance of rich clusters in 
the Universe at 2; = 0. The parameters controlling star formation and supernova feedback are 
set so that galaxies with circular velocities of 220 km s^^ have an I-band magnitude consistent 
with the Tully-Fisher relation of Giovanelli et al (1997). 

We show that different ways of treating star formation and feedback in the models have 
a strong influence both on the galaxy luminosity function and on the slope and amplitude 
of the two-point correlation function. In order for the high-density rCDM model to come 
close to matching the observations, we have to suppress the formation of galaxies in low-mass 
halos in the field. This is achieved by assuming that supernova feedback is so efficient that 
it can eject a substantial fraction of the baryons from the potential wells of low-mass dark 
matter halos. Even so, the model still produces an excess of very bright galaxies in groups 
and clusters. The low-density ACDM model produces too /eto star-forming field galaxies, even 
if feedback is weak and gas never escapes from dark halos. We have also investigated the 
effect of dust extinction on predictions for the colour distributions and clustering properties 
of galaxies. Given the uncertainties in modelling some of the critical physical processes, we 
conclude that we cannot reliably constrain the values of cosmological parameters using the 
properties of the galaxy distribution at 2; = 0. 
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2 The simulations 



The simulations we use were run as a part of the GIF project, a joint effort of astrophysicists 
from Germany and Israel. Its primary goal is to study the formation and evolution of galaxies 
in a cosmological context using semi-analytical galaxy formation models embedded in large 
high-resolution A^-body simulations. The code used for the GIF simulations is called Hydra. 
It is a parallel adaptive particle-particle particle-mesh (AP^M) code (for details on the code, 
see Couchman, Thomas, & Pearce 1995; Pearce & Couchman 1997). The current version was 
developed as part of the VIRGO supercomputing project and was kindly made available by 
them for the GIF project. The simulations were started on the CRAY T3D at the Computer 
Centre of the Max-Planck Society in Garching (RZG) on 128 processors. Once the clustering 
strength required a finer base mesh and so an even larger amount of total memory, they were 
transferred to the T3D at the Edinburgh Parallel Computer Centre (EPCC) and finished on 
256 processors. 

A set of four simulations with N — 256^ and with different cosmological parameters was 
run (Table 1). All the simulations are "cluster normalized". White, Efstathiou, & Frenk 
(1993) introduced this way of fixing the amplitude by determining ag, the square root of 
the variance of the density field smoothed over 8 h^^ Mpc spheres, such that the observed 
abundance of high-mass clusters is matched. More recent determinations of erg using the 
observed cluster X-ray temperature function (Eke, Cole, & Frenk 1996; Viana & Liddle 1996) 
yield similar results. For the low-density GIF simulations, the result by Eke et al. (1996) was 
taken. For the Q = 1 simulations, slightly larger values than suggested by Eke et al. (1996) 
were adopted, according to the earlier result by White et al. (1993). 

The parameters shown in Table 1 were chosen not only to fulfil cosmological constraints, 
but also to allow a detailed study of the clustering properties at very early redshifts. The 
masses of individual particles are 1.0 x 10^° Mq and 1.4 x 10^° Mq for the high- and low- 
Vt models, respectively. The gravitational softening was taken to be 30 kpc. Gravitational 
lensing by clusters in this set of four models has been studied by Bartelmann et al (1998). 

In this paper we focus on two of the models: rCDM and ACDM. Both have shape param- 
eter r = 0.21. A value of F = 0.21 is usually preferred by analyses of galaxy clustering, for 
example Peacock & Dodds (1994). This is achieved in the rCDM model, despite = 1 and 
h = 0.5, by assuming that a massive neutrino (usually taken to be the r neutrino) was present 
during the very early evolution of the Universe and came to dominate the energy density for 
a short period. It then decayed into lighter neutrinos which are still relativistic, thus delaying 
the epoch when matter again started to dominate over radiation. The neutrino mass and 
lifetime are chosen such that F = 0.21. For a detailed description of such a model see White, 
Gelmini, & Silk (1995). 

3 Construction of Halo Merger Trees and Identification 
of Galaxy Positions 

In order to follow the merging history of dark matter halos in the simulations, we store particle 
positions and velocities at 50 different output times, spaced in equal logarithmic intervals in 
redshift from z — 20 to z — 0. The construction of merger trees from these simulation outputs 



4 



involves the following steps. 

A friends-of-friends group-finding program is used to locate virialized halos. We adopt a 
linking length which is 0.2 times the mean interparticle separation. Only halos containing at 
least 10 particles are included in our halo "catalogues" . Tests show that 10-particle halos are 
stable systems. More than 95 % of 10 particle halos identified at one output time are still 
located within groups of 10 particles or more at subsequent times. By "located", we mean 
that more than 80% of the particles in one halo are present in the same halo at the later time. 
Halos with masses below about 7 particles do not survive over many output times according 
to this criterion. 

The lowest luminosity galaxy that we are able to resolve in these simulations thus corre- 
sponds to a galaxy in a halo with 10 particles or ~ 10^^h~^MQ. In these models, the Milky 
Way halo has a mass of ~ 2 x 10^^/i~^Mq, so the faintest resolved galaxies are roughly a tenth 
as bright (i.e. comparable in luminosity to the Large Magellanic Cloud). 

We then compute a set of physical quantities for all the halos in the catalogues. These 
are: 

1. the central particle index. This is the index of the most-bound particle in the halo. This 
particle has particular significance as it marks the position of the central galaxy of the 
halo, ie the galaxy onto which gas cools and where it forms stars. 

2. Rvir, the virial radius, defined as the distance from the central particle within which the 
overdensity of dark matter is 200 times the background density. 

3. Myir, the virial mass. This is the mass of dark matter contained within Ryir- 

4. Vc, the circular velocity {Vc — {GMyir/RvirY^'^). 

We begin with the first simulation output that has at least one halo with 10 particles or 
more. The central particles of these halos mark the locations of the first galaxies identified 
in the simulation. We then go to the next output time and loop through all the halos in the 
catalogue, searching for progenitor halos at the previous time. A halo at redshift Zi is defined 
to be a progenitor of a halo at redshift Zq < z, if 

1. more than half its particles are included in the halo at redshift zq] and 

2. its central particle is also included in that halo. 

The most massive progenitor of a halo has a special status - the properties of its central 
galaxy are transferred to the central galaxy of the new halo. We thus reposition the central 
galaxy at each output time. The particle index corresponding to the central galaxy changes, 
but its associated mass and luminosity evolve in a smooth fashion. This procedure ensures 
that cold gas always settles at the centre of a halo. The central galaxies of less massive 
progenitors become satellites. The particle index of a satellite galaxy then remains fixed until 
the present day. A satellite is said to "belong" to a halo if its particle index is among those 
linked together by the groupfinder. 

This procedure is repeated at every output time until z — 0. At each timestep, and for all 
the halos in the associated catalogue, we store the following information: 
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1. The index of the central particle of the largest progenitor of the halo (zero if the halo 
had no progenitors). 

2. The indices of all satellite galaxies contained within the halo (particles that were central 
galaxies of smaller halos at earlier times, and that are now incorporated within the 
present halo) 

This is sufficient to follow the evolution and merging of all galaxies within the simulation 
volume using the scheme outlined in the next section. 

Occasionally a galaxy identified at one time is not included within any halo at a subsequent 
time. This may happen for the following reasons: 

1. A halo that was above the 10-particle resolution limit at one time, may fall below this 
limit at a later time. In this case, there is a central galaxy that does not belong to any 
halo in the next simulation output. 

2. Satellite galaxies are occasionally "ejected" out of halos, particular during mergers. In 
most cases, the satellite will fall back into the halo at a later time. 

We keep an index list of these "lost" galaxies, checking at each subsequent output time to see 
whether they have been re-incorporated into a halo. If the recovered particle was previously 
a central galaxy, its properties are transferred to the central galaxy of the new halo if the 
difference in mass between the old and new halos is small (less than a factor 2). Otherwise, 
the recovered particle becomes a satellite galaxy within the new halo. Recovered satellite 
galaxies simply join the satellite population of the new halo. Note that lost galaxies comprise 
only a few percent of the total galaxy population at all times. Moreover, most of them are 
faint. 

4 The Physical Processes Governing Galaxy Formation 

Our treatment of the physical processes governing galaxy formation is very similar to that 
described in Kaufi'mann, White & Guiderdoni (1993, KWG) and in White & Frcnk (1991). 
The reader is referred to these papers for more detailed discussion and for derivations of some 
of the equations. Some of our prescriptions, in particular those for feedback and for galaxy- 
galaxy merging, have changed since these papers were published. We have also incorporated a 
number of extra features, such as dust extinction and starbursts during galaxy-galaxy mergers. 
For completeness, we now present a summary of our prescriptions for the various physical 
processes. A detailed discussion of how these are implemented is given in section 4.8. 

4.1 Gas Cooling 

We adopt the simple model for cooling first introduced by White & Frenk (1991). For sim- 
plicity, dark halos are modelled as isothermal spheres truncated at their virial radius Rmr- 
We assume that the hot gas always has a distribution that exactly parallels that of the dark 
matter. The total mass of hot gas in the halos is given by 

Mhot = ^bMyir - M, - M^ool - Meject: (1) 
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where is the baryon density of the Universe, Myir is the virial mass of the halo, is the 
total mass of stars that have formed in the halo, Mcooi is the total cold gas contained in the 
halo and Meject is the mass of gas ejected out of the halo by supernova feedback (see section 
4.3). 

The gas temperature can then be derived from the circular velocity of the halo using the 
equation of hydrostatic equilibrium: 

T = 35.9(K/kms-^)2K (2) 

At each radius in the halo, we define a local cooling time through the ratio of the specific 
energy content to the cooling rate A{t) . 

tcoolv) = 7: 2/ \ \ u\ \ ) 

where Pg{r) is the gas density, ne(r) is the electron density, nip is the proton mass and /imp is 
the molecular weight of the gas. In this paper, we do not include chemical enrichment. The 
gas is assumed to have solar metallicity at all times and we use the solar metallicity cooling 
curve in figure 9.9 of Binney & Tremaine (1987). (see Kauffmann 1996b and Kauffmann & 
Chariot 1998 for a description of models including chemical evolution). 

At a given redshift z, a cooling radius rcooi can be defined as the radius within the halo 
where the cooling time is equal to the age of the Universe. For the case of an Einstein-de 
Sitter cosmology, 

tcooi{rcooi)^lH^\l + z)-'/'. (4) 

At high redshifts, for small halos, and for high gas fractions, the coohng radius is larger 
than the virial radius of the halo. We assume that in the absence of supernovae, all the hot 
gas in the halo would settle to the centre on a timescale given by the halo dynamical time 
{Ryir/Vc). We thus write the instantaneous cooling rate as 

M^oi{Vc,z) = ^%^. (5) 

. The gas content of the halo is further affected by feedback (see next section) and by the 
infall of new material, which we determine directly from the merging trees. 

At later times, for larger halos, and for low gas fractions, the cooling radius lies inside the 
virial radius and the rate at which gas cools is calculated using the equation 

Mcooi(Vc,z) ^^T^Pg{rcooi)rloor-^- (6) 

As has been shown in KWG and in Aragon-Salamanca, Baugh & Kauffmann(1998), for 
massive halos the cooling rates given by equation 6 lead to the formation of central cluster 
galaxies that are too bright and too blue to be consistent with observation if the cooling gas is 
assumed to form stars with a standard initial mass function. It should be noted that cooling 
flows of hundreds of solar masses per year are observed in a number of clusters (see for example 
Fabian, Nulsen & Canizares 1991; Allen & Fabian 1997), but that the fate of the cooling gas 
remains a mystery. One hypothesis is that it may condense into cold clouds instead of stars 



7 



(Ferland, Fabian & Johnstone 1994). In our models, we assume the gas coohng in halos with 
Vc > 350 km s~^ docs not form visible stars. Note that this is the lowest circular velocity 
at which we can suppress star formation in cooling flows without destroying our fit to the 
TuUy-Fisher relation (see section 6.2). 

4.2 Star Formation 

As in KWG, we adopt a simple star formation law of the form 

= aM^oid/tdyn: (7) 

where a is a free parameter and t^j/n is the dynamical time of the galaxy. For a central galaxy 
in a halo, the dynamical time is given by 

tdyn = O-^Rvir/Vc- (8) 

For disk galaxies, this is motivated by noting that, if gas collapses to a centrifugally supported 
state within an isothermal halo while conserving angular momentum, the contraction factor 
is ~ 2A, where A is the spin parameter of the gas, assumed to be the same as for the halo. 
N-body simulations find that A scatters around a value of ~ 0.05. For a satellite galaxy, tdyn 
is held fixed at the value when the galaxy was last a central galaxy. 

It should be noted that according to the simple spherical collapse model, the virial radius 
of a dark matter halo scales with circular velocity and with redshift as R^ir oc Vc{l + z)~^^^, 
reflecting the fact that halos are smaller and denser at earlier epochs. This means that tdyn 
is independent of the circular velocity of the halo, but will decrease at higher redshift, so star 
formation rates are higher in halos of the same cold gas content at high z. 

The star formation law in equation 7 has received considerable empirical support from 
a recent study of the star formation rates and gas masses in 61 nearby spiral galaxies and 
36 "starburst" systems by Kennicutt (1997), who finds that such a law can fit the data over 
several orders of magnitude in star formation rate and gas density. In normal spirals, about 
10 percent of the available gas is turned into stars per orbital time. 

4.3 Feedback from Supernovae 

The effect of energy ejected by supernova explosions into the interstellar medium of a galaxy 
has profound implications for the observed properties of galaxies. As shown by Cole et al 
(1994) and Somerville & Primack (1998), strong feedback in low-mass galaxies is required 
to fit the fiat {a > —1.3) faint-end slope of the galaxy luminosity function. Kauffmann & 
Chariot (1998) have demonstrated that substantial feedback in massive (~ L*) galaxies is 
also needed to fit the observed slope of the colour-magnitude relation of elliptical galaxies. 
Unfortunately, both theoretical and observational understanding of how feedback operates in 
different types of galaxies is extremely limited at present, so we have no option other than to 
experiment with a variety of different prescriptions in order to see what difference they make 
to our results. 

Using basic energy-conservation arguments, it is possible to estimate how much cold gas 
could be reheated to the virial temperature of the halo for a given mass of stars formed in 
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the galaxy. For the Scalo (1986) initial mass function we use to model the evolution of the 
stellar populations in our galaxies (see section 4.6), the number of supernovae expected per 
solar mass of stars formed is tjsn = 5 x 10~^Mq^. The kinetic energy of the ejecta from each 
supernova, Esn, is about 10^^ erg. If a fraction e of this energy is used to reheat cold gas to 
the virial temperature of the halo, the amount of cold gas lost in time At can be estimated as 

AM..,..,^4 '^-"f.f^" At. (9) 
3 

In the models, e is treated free parameter. 

One major uncertainty is whether the gas heated by supernova explosions will leave the 
halo. In our previous work, reheated gas was always retained within the halo. Another 
possibihty is that reheated gas will be ejected out of halo. In the models of Cole et al (1994), 
reheated gas was removed until the halo grew in mass by a factor of two or more, whereupon 
it was added once again to the hot gas component. The re-incorporation of expelled gas at a 
later time ensured that the total baryonic mass in halos was conserved and that the baryon 
fraction in clusters was close to the global value. 

In this paper, we will experiment with both feedback prescriptions. We call the model in 
which reheated gas is always trapped within the halo the "retention" model. The model in 
which reheated gas is expelled from the halo will be called the "ejection" model. According to 
our star formation prescription, stars form efficiently in low mass halos at high redshifts. In 
the ejection model, the energy injected into the ISM by the first generation of star formation 
is sufficient to expel most of the gas from the halo. The star formation rate then drops and the 
galaxy fades until the gas is re-incorporated on the next collapse. Stars in low-mass galaxies 
are thus formed in a series of "bursts" associated with each factor 2 doubling in halo mass. 
This bursting behaviour is much less pronounced for massive galaxies, since their potential 
wells are deeper and much less gas will be expelled (equation 9). 

The effect of the two feedback prescriptions on the typical star formation history of a 
galaxy (more specifically, the central galaxy of a halo with Vc — 220 km s~^) is shown in figure 
1. As can be seen, the ejection prescription results in both higher and more irregular rates of 
star formation at high redshift when the galaxies reside within low-Vc halos. At low redshifts, 
the halo potential wells are deeper and so more effective at retaining the gas heated during 
supernova explosions. The star formation rates for for the two prescriptions thus do not differ 
very much at the present epoch. 



4.4 Merging of Galaxies 

N-body plus smoothed particle hydrodynamic (SPH) simulations of the assembly of galaxies 
in a hierarchical Universe show that as dark matter halos coalesce, the embedded disk galaxies 
merge on a timescalc that is consistent with dynamical friction estimates based on their total 
(gas + surrounding dark matter) mass (Navarro, Frenk & White 1995). The dynamical friction 
timescale, in the form given by Binney & Tremaine (1987) is: 
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Figure 1: A comparison of the star formation histories for 3 central galaxies in halos of K = 
220 kms~^ for the retention and the ejection feedback schemes. 

where Vc is the circular velocity of a singular isothermal sphere representing the primary 
halo, Mgat is the mass of the orbiting satellite, ln(A) is the Coulomb logarithm which we 
approximate as In^Myir/Mgat)- The function /(e) allows for the angular momentum of the 
satellite's orbit, expressed in terms of the circularity parameter e = J/ Jc{E), the ratio of the 
angular momentum of the satellite to that of a circular orbit with the same energy. Lacey & 
Cole (1993) show that /(e) ~ e"'^^ for e > 0.02. Finally, rc{E) is the radius of the circular 
orbit with the same energy as the satellite orbit. Navarro, Frcnk & White show that equation 
10 provides, on average, a fair estimate of the merger timescalc of a satellite galaxy, provided 
Mgat is taken to be the total gaseous + dark matter mass of the satellite. 

In our models, we use equation 10 to compute the timescale for an accreted satellite to 
reach the centre of the halo and merge with the central galaxy. The orbital eccentricity 
parameter e is drawn randomly from a uniform distribution from to 1. The radius Vc is set 
equal to Rvir, the virial radius of the primary halo. Mgat is taken as the baryonic mass of the 
satellite plus the mass of its surrounding halo at the time it was last a central galaxy. 

Note that it is assumed that satellite galaxies merge only with the central galaxy of the 
primary halo and only after time Tdynf- If the primary halo is later accreted by a larger 
system, new orbital parameters for all the remaining unmerged satellites are drawn, the dy- 
namical friction timescales are recomputed and the merger clock is reset to zero. In practice, 
substructure in halos will not be erased immediately and for a while a satellite galaxy may 
still be able to merge with its old central object. 
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4.5 Formation of Elliptical Galaxies and Spiral Bulges 

If two galaxies merge and the mass ratio between the sateUite and the central object is 
greater than 0.3, we add the stars of both objects together and create a bulge component. If 
Mgat/Mcentrai < 0.3, wc add the stars and cold gas of the satellite to the disk component of 
the central galaxy. The value of 0.3 is motivated by a series N-body simulations of merging 
disk galaxies of unequal mass (Barnes, private communication). When a bulge is formed by 
a merger, all cold gas present in the two galaxies is transformed into stars in a "starburst" 
with a timescale of 10^ years. Further coohng of gas in the halo may lead to the formation of 
a new disk. 

The morphological classification of galaxies is made according their B-band disk-to-bulge 
ratios (Simien & de Vaucouleurs 1986). If M{B)huige — M{B)totai < 1 mag, then the galaxy is 
classified as early- type (eUiptical or SO). 

It should be noted that although we track the formation of galaxies in halos as small 
as 10 particles, we do not accurately predict the morphologies of galaxies contained in such 
halos, simply bccaiisc their merging histories are not resolved. Accurate morphologies are only 
obtained for central galaxies in halos with ~ 100 particles, i.e. for galaxies with luminosities 

4.6 Stellar Population Models 

We use the new stellar population synthesis models of Bruzual & Chariot (in preparation), 
which include updated stellar evolutionary tracks and new spectral libraries. In this paper, 
all stars are assumed to have solar metallicity. The star formation history of any galaxy 
can be approximated by a series of delta-function "bursts" of different masses. The stellar 
population models are used to generate lookup tables of the luminosity of a burst of fixed 
mass as a function of age in each photometric band. The magnitude of the galaxy at 2; = is 
calculated by summing the mass- weighted luminosities of each burst. 

We have adopted a Scalo (1986) initial mass function with upper and lower mass cutoffs 
of 100 Mq and 0.1 Mq. As shown in Kauffmann & Chariot (1998), the stellar mass-to-light 
ratios of an old (8-10 Gyr) stellar population are then in good agreement with the observed 
mass-to-light ratios of elliptical galaxies within an effective radius. 

4.7 Dust Extinction 

In the models there is a strong distinction between central and satellite galaxies. Gas cools 
continuously onto central galaxies and they form stars at a roughly constant rate. Satel- 
lite galaxies lose their supply of new gas and, as a result, have exponentially declining star 
formation rates, leading to redder colours and lower gas fractions than central galaxies. A 
central galaxy is thus considerably brighter than a satellite companion of the same stellar 
mass at ultraviolet and optical wavelengths. As we will show in section 6.4, this infiuences 
the amplitude of the correlation function of galaxies selected according to B-band magnitude. 

It is a well-known but oft-ignored observational fact that gas-rich, star-forming galaxies 
contain dust, which absorbs a substantial fraction of the light emitted in the short wavelength 
part of the spectrum and re-radiates it in the far-infrared. We include a simple, empirically- 
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motivated recipe for dust extinction in our models in order to investigate how much this effect 
can influence our estimates of clustering in the models. 

Wang & Heckman (1996) have studied the correlation of the optical depth of dust in 
galactic disks with the total luminosity of the galaxy. They studied 150 normal late-type 
galaxies with measured far-ultraviolet (UV, A ~ 2000A) fluxes and compiled the corresponding 
far-infrared (FIR, A ~ 40 — 120/im) fluxes measured by the IRAS satellite. They then modelled 
the absorption and emission of radiation by dust with a simple model of a uniform plane- 
parallel slab in which the dust that radiates in the IRAS band is heated exclusively by UV 
hght from nearby hot stars. They find that their observed UV-to-FIR ratios can be explained 
by the face-on extinction optical depth r varying with the intrinsic UV luminosity of the 
galaxy as 

T OC To{L/L,f, (11) 

with /3 ~ 0.5 The same scaling law was also able to account for the Hf3/Ha ratios measured 
for a subset of the galaxies. Re-expressed in the blue-band, Wang & Heckman derive a total 
extinction optical depth of tb,* = 0.8 ± 0.3 at the fiducial observed blue luminosity of a 
Schechter galaxy (M,(5) ='-19.6 + 5 log /i). 

We use the Galactic extinction curve of Cardelli, Clayton & Mathis (1989) to derive t\/tb- 
Although the extinction curves of smaller, less-metal rich galaxies such as the LMC and SMC 
are quite different at ultraviolet wavelengths, at optical wavelengths they are all fairly similar. 
This leads to face-on extinction for an disk galaxy of 0.9, 0.8,0.4 and 0.08 mag in the B,V,I 
and K bands respectively. 

Finally, one must also take into account the inclination of the galaxy to the hne-of-sight 
when making an extinction correction. For a thin disk where dust and stars are uniformly 
mixed, the total extinction in magnitudes is 

( \ — p-r\secO\ 

AA = -2.51og — . (12) 

\ Tx sec 9 J 

In our dust-corrected models, we apply an extinction correction to every galaxy forming stars 
at a rate greater than 0.5 Moyr"^. We use M^(B) — —19.6 -|- 5 log h — 0.9 in equation 11, and 
scale Tx using the uncorrected B-magnitudes predicted by the model. We then pick a random 
inclination for each galaxy and apply the final correction using equation 12. 

4.8 Detailed Implementation of the Prescriptions 

In the simulation, each galaxy carries a number of "labels" corresponding to different physical 
properties. These are: 

1. M*, total stellar mass 

2. M^{bulge), stellar mass of the bulge component 

3. Mcooh cold gas mass 

4. L{B, V, R, I, K...), total present-day luminosity in a given waveband 

5. Lhuige{B, V, R, I, K...), the present-day bulge luminosity 
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Satellite galaxies carry two additional labels: 

1- truer g-i the time until the satellite galaxy should merge with the central object. 
2. «merg, the iudcx of the halo in which the satellite resided at the previous timestep. 

In the ejection feedback scheme, each halo carries an array Meject{Mprog) , which is the 
amount of gas ejected by galaxies in progenitor halos of mass Mp^og- 

At each output time, we loop over all the halos in the catalogue and follow a series of 
steps: 

1. We re-incorporate gas that was ejected out of progenitor halos more than a factor of 
two less massive than the present halo 

2. We calculate the mass of hot gas in the halo at the start of the timestep. This is given 

by 

Mhot = ^bMyir-Y,{Mcool + M,)- MejectiKrog), (13) 

gals Mprog 

where J2gais is the sum over all galaxies in the halo. Note that this assumes that when 
halos merge with each other, all gas that is not already cooled is shock heated to the virial 
temperature of the new halo. The sum over Mp^og extends only over those progenitors 
more massive than Myir/2. 

3. We loop through the satellite galaxies in the halo and update their merging timcscales. 
If the satellite was not in the largest progenitor of the halo at the previous timestep, 
its merging clock is reset. If the galaxy is predicted to merge during the timestep, it 
is flagged. The time between the beginning of the timestep and the merging event is 
noted. 

4. We then solve a set of coupled differential equations for the time evolution of the cold 
gas and stars in each of the galaxies in the halo. We adopt a small timestep for this 
calculation (there are typically 100 timesteps between each pair of simulation outputs). 
For central galaxies, the change of the cold gas component over time At is given by: 

Mcoolit + At) = Mcoolit) - (M* + Mreheat - Mcool)At, (14) 

where Mcooh ^* ^'^^ ^reheat ^rc given by equations 5,6 and 7(or 9) respectively. For 
satellite galaxies, the equation is the same except there is no source term from gas cooling 
in the halo. Infall of cold gas onto satellites is assumed to be disrupted as soon as they 
are accreted so that star formation within them continues only until their existing cold 
gas reservoirs are exhausted. 

The stars formed during time At are added to M*. We look up the present-day lumi- 
nosity of a burst of mass M^At and age t in our population synthesis tables and update 
the luminosities of the galaxies. Finally, the hot gas mass of the halo at the end of the 
timestep is calculated using 

Mhot{t + At) = Mhot{t) - Yi^cool - Mreheat) At. (15) 

gals 

Note that this is valid for retention feedback. For ejection feedback, reheated gas is 
added instead to Meject(-^wr); where M^j^ is the virial mass of the halo. 
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5. If a satellite is predicted to merge during the timestep At, its stars and gas are added 
to the central galaxy. If the mass of the satellite is greater than one-third the mass of 
the central object, we set the stellar mass of the bulge component equal to the mass of 
the central galaxy plus the mass of the merged satellite. All the cold gas in the merger 
remnant is converted into stars at a constant rate over 10^ years. Stars that form in 
such a burst do not reheat or eject cold gas. In addition, we assume that the the burst 
is decoupled from the evolution of rest of the galaxy, i.e. we set Mcooi = in the above 
equations immediately after the merger takes place. Both assumptions are motivated by 
N-body simulations, which show that gas in merging disk galaxies quickly loses angular 
momentum and ends up in a dense knot at the very centre of the remnant, and that 
feedback has rather little effect (Barnes & Hernquist 1996). 

4.9 Normalizing the Models 

For a given set of cosmological initial conditions (including the baryon density Qf,), there 
are two "free parameters" in our models - the star formation efficiency a and the feedback 
efficiency e. As in KWG, we tune these two parameters to match the luminosity and cold 
gas mass of a fiducial reference galaxy, which we take to be the central galaxy in a halo with 
circular velocity Vc — 220 km s~^. In most previous work we have normalized the models using 
standard parameters for the Milky Way, but here, following the suggestions of Somerville & 
Primack (1998), we adopt a normalization that is consistent with the velocity-based zero-point 
of the I-band Tully-Fisher relation as determined by Giovanelli et al (1997). Matching to I- 
band rather than B-band data ought to be more robust, since the I-band is considerably less 
affected by dust extinction and star formation. Moreover, a very large amount of data on the 
I-band Tully-Fisher relation is now available. Giovanelli et al obtain the following fit, based 
on 555 galaxies in 24 clusters: 

Mj-51og/i = -21.00 ±0.02 - 7.68 ±0.14(logiy- 2.5) (16) 

To convert between the measured HI line-widths W and the model circular velocities 
we assume W = 2Vc. In practice, there is substantial uncertainty in the relation between 
the circular velocity of the halo and that of the disk (e.g. Mo, Mao & White 1998), since 
the transformation between the two quantities depends on the detailed density profile of the 
halo and whether or not gas loses angular momentum before settling onto the disk. We 
set the parameters a and e so that the central galaxy in a halo with Vc = 220 km s~^ 
has Mj — 51og/i ~ —22.1 and a cold gas mass ~ IO^^Mq. It should be noted that this 
normalization is considerably brighter than that adopted by KWG, who took the observed 
B-band magnitude of the Milky Way to be 5 ~ —20.5. For a Hubble constant of 50 km s~^ 
Mpc~^ and assuming S — / ~ 1.8 for a spiral galaxy (De Jong 1996), the new normalization 
means that the central galaxy in a halo with Vc = 220 km s~^ has a B-band magnitude of 
-21.8 - more than a factor 3 more luminous than the observed value. Gonsistency of this 
Tully-Fisher normalization with the estimated magnitudes of our own Galaxy and of M31 
is only obtained for substantially higher values of the Hubble constant, which do not give 
acceptable ages for a Universe with Q — 1. 

The fact that there are only two free parameters in the model should not be interpreted 
as an indication that the physical processes governing galaxy formation are well-specified. 
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On the contrary, we will show that the way in which we choose to implement feedback and 
whether or not dust is included in the model can have a very large effect both on the luminosity 
function and on the amplitude of the correlation function on scales below a few Mpc. We 
will attempt to clarify how different parametrizations of these processes affect our results and 
which processes appear to be most critical if one is to succeed in matching the observational 
data. 

5 Comparison of galaxy properties for semi-analytic and 
simulation merger trees 

If halo merging trees derived using methods based on the extended Press-Schechter formalism 
(Cole 1991; Kauffmann & White 1993; Rodrigues & Thomas 1996 Somerville & Kolatt 1998) 
were statistically equivalent to the trees found in the N-body simulations, the same galaxy 
luminosity functions would be obtained using the two approaches, provided the same recipes 
for coohng, star formation and feedback were employed, and the same halo mass resolution 
limit was adopted. In practice, the halo mass function derived by the Press-Schechter argu- 
ment does not fit very accurately the mass function in N-body simulations. For halo masses 
between 10^^ and 10^"^ M©, the Press-Schechter theory predicts roughly twice as many halos 
as are actually found in the GIF simulations for both the rCDM and the ACDM cosmologies. 
Most of the missing mass in the simulations is in the form of "unresolved" material - single 
particles or groups with less than 10 members. This discrepancy between the Press-Schechter 
and simulation mass functions propagates to the halo progenitor distributions and to the 
Monte-Carlo merging trees (Somerville, Lemson & Kolatt, in preparation). 

Here we compare the the predictions of the two approaches for the luminosity functions 
of galaxies within halos of given mass and for the "field" luminosity function. In figure 2, we 
plot the number of galaxies brighter than a given V-magnitude in halos with circular velocity 
Vc- Solid squares show mean results derived using analytic trees constructed as described 
in Kauffmann & White (1993). Solid circles show results derived from the simulation. The 
error bars on the points show the rms scatter in the number of galaxies in halos of given 
circular velocity. We perform the comparison using the rCDM model and we use the same 
prescriptions (from KWG) for cooling, star formation, feedback, merging and stellar popula- 
tion synthesis in both cases. We truncate the analytic trees at a halo mass of 2 x lO^^M© so 
that the resolution is the same as in the simulations. The four panels in figure 2 show the 
number of galaxies per halo brighter than My — —18, —19, —20 and —21. As can be seen, the 
mean number of lower-luminosity galaxies per halo (My > —20) agrees remarkably well. The 
analytic approach appears to overpredict the number of bright galaxies in high halos by a 
factor of ~ 2, but since there are not many such halos in the simulation and the number of 
bright galaxies within them is small, the discrepancy may not be statistically significant. 

Figure 3 compares the field luminosity functions derived using the semi-analytic trees 
and the simulations. The top panel shows the difference between the abundance of halos 
predicted by Press-Schechter theory, and that found in the simulation. The Press-Schechter 
overprediction of halo abundance at intermediate masses in quite evident. The middle panel 
compares the magnitude of the central galaxy as a function of halo circular velocity. There 
is good agreement between the two approaches. Note that the dip in magnitude at Vc — 500 
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Figure 2: The number of galaxies brighter than a given V-magnitude in halos with circular velocity 
Vc- Solid circles show results derived using the analytic trees constructed as described in KaufFmann 
&: White (1993). Solid squares show results derived from the simulation. Error bars show the rms 
deviation between the number of galaxies in different halos. 



16 



km arises because we assume stars do not form in the cooling flows of halos with circular 
velocities greater than this value. The bottom panel compares the V-band "field" luminosity 
functions. Except perhaps near the "knee", the agreement between the two approaches is 
excellent. 

The comparison in this section demonstrates that although merging trees derived using 
the extended Press-Schechter theory differ in detail from those found in the simulations, the 
two methods give very similar results both for the mean number and luminosity of galaxies 
within halos, and for the scatter in these quantities. 

6 Results of the Models 

6.1 Slices from the simulation 

In figure 4, we show slices of thickness 8 Mpc from the rCDM simulation. The distribution 
of dark matter is shown in the top left panel. In the top right panel, all the galaxies in the 
slice with magnitudes brighter than M{B) = — 19.0 + 5 log /i are plotted as sohd white circles. 
Galaxies with magnitudes brighter than M{B) — —17.5 + 51og/i are shown in the bottom 
panels. We divide these galaxies into two equal subsamples at the median star formation rate 
per unit stellar mass. The bottom left panel shows the distribution of galaxies with low star 
formation rates. The bottom right panel shows galaxies with high star formation rates. 

Galaxies of all luminosities trace the filaments and knots visible in the dark matter dis- 
tribution. Galaxies with high star formation rates are less clustered than galaxies with low 
star formation rates. Star-forming galaxies tend to avoid dense clusters and groups and occur 
more frequently in voids. 

Figure 5 shows a similar set of slices from the ACDM simulation. The area and thickness 
of these slices have been chosen to have the same dimensions in redshift space as those in figure 
4 (85/i~^ Mpc x85h~^ Mpc x8h~^ Mpc). The same systematic effects are visible in the two 
figures and the most obvious difference between them is the difference in galaxy abundance. 
We return to this in section 6.3. 

6.2 The Tully-Fisher Relation 

As described in section 4.9, all models are normalized so that the average I-band magnitude 
of the central galaxy in a halo of circular velocity Vc = 220 km s^^ is —22.1 + 51og/i, in 
accordance with the zero point of the I-band Tully-Fisher relation derived by Giovanelli et 
al. (1997). Figure 6 shows the Tully-Fisher relation for spiral galaxies in the simulations. 
We have selected central galaxies with 1.5 < M{B)fyuige — M{B)tot < 2.2 (appropriate for 
Sb/Sc type galaxies). The solid line is the fit to the data (equation 16). The upper two 
panels show the rCDM model for the two different feedback prescriptions discussed in section 
4.3. Although the slopes are similar, the ejection prescription results in considerably more 
scatter, particularly at low velocity widths, where there is a marked tail of low-luminosity 
galaxies. The lower panel shows the ACDM model with retention feedback. We do not show 
a ACDM model with ejection feedback as this model would result in too low a total luminosity 
density, as explained in the next section. It is interesting that the scatter obtained for both 
the ACDM and rCDM relations, even for retention feedback, is not much smaller than the 
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Figure 3: Top panel: A comparison of the abundance of iialos as a function of circular velocity 
predicted by the Press-Schechter theory (solid) and the abundance derived from the simulation 
(dotted). Middle panel: Magnitude of the central galaxy as a function of halo circular velocity. 
Error bars show the rms scatter between halos. The Prcss-Schcchter theory halos are solid and 
the simulation is dotted. Bottom panel: The V-band field luminosity functions derived using the 
Press-Schechter approach (solid) and the simulations (dotted). 

Figure 4: Slices of thickness 8 h^^ Mpc from the rCDM simulation. The distribution of 
dark matter is shown in the top left panel. In the top right panel, all galaxies with M{B) < 
—19.0 + 51og/i are plotted as solid white circles. Galaxies with M{B) < —17.5 + 51og/i 
are shown in the bottom panels. We divide these galaxies into two equal subsamplcs at the 
median star formation rate per unit stellar mass. The bottom left panel shows the distribution 
of galaxies with low star formation rates. The bottom right panel shows galaxies with high 
star formation rates. 

Figure 5: Slices of thickness 8 Mpc and length 85 Mpc from the ACDM simulation. 
The panels are as described in figure 4. 
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Figure 6: The I-band TuUy-Fisher relation for spiral galaxies in the simulations. The solid lines 
show the relation of Giovanelli et al (1997). 

observed scatter. The model scatter arises because of differences in star formation history 
and stellar mass between central galaxies in different halos of the same circular velocity. As 
pointed out by Mo, Mao & White (1998), scatter in the spin parameter of dark matter halos 
leads to a spread in contraction factors for disks forming in halos of a given circular velocity. 
Halos with high values of A produce disks with rotation curves that are slowly rising, whereas 
halos with low values of A produce disks with rotation curves that rise steeply and then decline. 
This leads to intrinsic scatter in the predicted Tully-Fishcr relation even if disks in halos of 
given circular velocity are assumed to have fixed mass and mass-to-light ratio. This must 
be added to the scatter already visible in figure 6. More complex modelling is necessary to 
investigate these issues in detail. 

6.3 The Luminosity Function 

In figure 7, we show B-band luminosity functions for the rCDM simulation. A sequence of 
different models is shown in order to illustrate the effect of including different prescriptions for 
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feedback and dust extinction. The results are compared with the B-band luminosity functions 
derived from several recent redshift surveys: the APM-Stromlo survey (Loveday et al 1992), 
the Las Campanas Redshift Survey (LCRS) (Lin et al 1997) and the ESO-Slice survey (Zucca 
et al 1997). Note that we only plot the luminosity functions down to magnitudes where 
our simulations are "complete" (i.e. the luminosity of the central galaxy in a halo with 10 
particles). 

In the top left panel, we plot the luminosity function for a model with retention feedback 
and no dust extinction. In this case, the model luminosity function exceeds the observed one 
by more than a factor of 10 at all luminosities. The inclusion of dust brings the model more in 
line with the observations at luminosities around L*, but there are still far too many galaxies 
at both fainter and brighter magnitudes. The model with ejection feedback agrees much 
better with the observations at L^, and below. This is because ejection feedback suppresses 
the formation of bright galaxies in low mass halos. The first generation of stars formed when 
the halo collapses generates enough energy to empty the halo of gas and halt further star 
formation. This process is only efficient in low-circular velocity systems, so the excess of 
galaxies brighter than still remains a problem. Including dust in the ejection model (lower 
right panel) helps somewhat, but the predicted abundance of bright galaxies is still more than 
a factor 10 above the Schechter function representation of the observations. 

In figure 8, we compare the model K-band luminosity functions with the Schechter fits 
derived by Gardner et al (1997) and Szokoly et al (1998) for two independent surveys. Unhke 
B-band magnitudes, the K-band luminosities of galaxies depend on their total stellar masses, 
rather than their present-day star formation rates. Moreover, K-magnitudes are only weakly 
affected by dust extinction. The fact that the model K-band luminosity functions also dis- 
agree with the data at the bright end shows that there are too many massive galaxies in the 
simulations. 

B-band and K-band luminosity functions for the ACDM model are shown in figure 9. 
Retention feedback is assumed. Even so, this model produces a factor ~ 2 — 3 too few 
galaxies at magnitudes around L*. There is again an excess of very bright galaxies, but the 
problem is less severe than in rCDM. The same trends are apparent in the K-band. 

In summary, the luminosity functions for both models have the wrong shape, with or 
without dust extinction. Instead of declining exponentially at bright magnitudes, they exhibit 
a gentler turn down. In previous work, KWG and Cole et al (1994) obtained much better fits 
to the bright end of the luminosity function (see also figure 3). As explained in section 4.9, 
this is because they adopted a much fainter normalization. If one is to obtain an exponential 
cutoff at high luminosities, bright galaxies must reside primarily in very massive halos on the 
exponentially declining part of the mass function. The shape problem is considerably worse 
for the rCDM model than for the ACDM model. Although both models are normalized to 
reproduce the observed abundance of rich clusters, the rCDM model has a higher abundance 
of intermediate mass halos than the ACDM model. It is the luminous galaxies in these halos 
that create the excess at magnitudes brighter than L*. The rCDM model produces too many 
galaxies below unless gas is efficiently ejected out of low mass halos. The ACDM model, 
on the other hand, produces too few galaxies even if feedback is inefficient and if gas never 
escapes these halos. 
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Figure 7: The B-band luminosity function for galaxies in the rCDM simulation. The four panels 
show models with different assumptions about feedback and dust extinction. The simulation results 
are shown as solid squares. The lines are Schechter fits to B-band luminosity functions from recent 
redshift surveys: 1) ESO-Slice (dotted), 2) APM-Stromlo (short-dashed), 3) LCRS (long-dashed). 
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Figure 8: The K-band luminosity function for galaxies in the rCDM simulation. The two panels 
show models with different assumptions about feedback. Dust extinction is negligible in the K-band. 
The simulation results are shown as solid squares. The lines are Schechter fits to K-band luminosity 
functions derived by Gardner et al (dotted) and Szokoly et al (dashed) 
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Figure 9: The B and K-band luminosity functions for galaxies in the ACDM simulation. The lines 
are as described in the previous two figures. 
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6.4 The Two-point Correlation Function 



In figure 10, we show the galaxy two-point correlation function ^(r) for the rCDM simulation. 

We select galaxies with B-band luminosities M{B) < —20.5. A sequence of models is again 
shown to illustrate different prescriptions for feedback and dust extinction. The results are 
compared with the real space correlation function measured from the APM galaxy survey 
(Baugh 1996). 

For the model with retention feedback, ^(r) follows a power law of slope 7 ~ — 1.7 down to 
scales ~ 1 Mpc, below which it flattens and then turns over. This is in marked contrast to the 
observed ^{r), which follows the same power law down to scales below 100 kpc. The inclusion 
of dust extinction lowers the luminosity of bright star-forming field galaxies and raises the 
amplitude of the correlation function on small scales, but ^(r) still flattens and is a factor 
6 below the observed value at r ~ 100 kpc. As described in the previous section, ejection 
feedback suppresses bright star-forming galaxies in low circular velocity halos and brings the 
faint end of the luminosity function into better agreement with the observations. As seen in 
the bottom left panel of figure 10, ejection feedback also prevents ^(r) from fiattening below 
1 Mpc. The results shown in the bottom right panel of figure 10, for the model with ejection 
feedback and dust extinction, agree rather well with the observations. 

Figure 11 shows the galaxy two-point correlation function for the ACDM simulation. Re- 
tention feedback is assumed, ^(r) does not turn over on small scales and the model with dust 
is actually too steep compared to observations. 

What determines the slope and amplitude of the galaxy correlation functions in our mod- 
els? In a recent paper. Mo & White (1996) developed a simple analytic model for the grav- 
itational clustering of dark matter halos. The positions and formation times of halos were 
determined from the statistics of the initial linear density field and modifications caused by 
gravitationally-induced motions were treated using a spherical collapse approximation. Mo & 
White showed that on large scales, the halo autocorrelation function is simply proportional to 
the dark matter correlation function, with the constant of proportionality dependent on the 
masses of the selected halos. Because halos are spatially exclusive and have a finite radius, the 
halo-halo correlation function drops below the mass correlation function on scales comparable 
to the diameter of typical halos in the sample. From these results, it follows that if there is 
only one bright galaxy per halo, the galaxy correlation function will simply follow the form 
of the halo correlation function, with a flattening and turnover on scales less than ~ 1 Mpc. 
In order for the galaxy correlation to continue as a power law to small scales, a substantial 
fraction of bright galaxies must exist as groups within the same halo. 

This is illustrated quantitatively in flgure 12, where we plot the fraction of galaxies in the 
simulation with M[B) < —20.5 that occur in halos of a given circular velocity. The rCDM 
model with retention feedback has a much higher fraction of galaxies in low Vc halos than 
the rCDM model with ejection feedback. Low Vc halos typically contain only one bright 
galaxy (see flgure 2), whereas high Vc halos contain many bright galaxies that contribute to 
clustering amplitude on sub-megaparsec scales. Figure 12 shows that there is a relatively small 
fraction of bright galaxies in low Vc halos in the ACDM model, even with retention feedback. 
Recall that both rCDM and ACDM are normahzed to flt the observed number density of 
rich clusters at the present day. Because the mass density in the ACDM model is lower, the 
number density of low mass halos is also smaller. ACDM thus "naturally" has more galaxies 
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Figure 10: The two-point correlation function ^(r) for galaxies in the rCDM simulation. Solid 
lines show the simulation results for different prescriptions for feedback and dust extinction. The 
points with error bars are taken from Baugh (1995). 
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Figure 11: The two-point correlation function ^(r) for galaxies in the ACDM simulation (solid 
line). Points are observational data from the APM survey. 
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Figure 12: The fraction of galaxies in the simulation with M{B) < —20.5 that occur in halos of 
a given circular velocity. Upper panel: rCDM, retention feedback (solid), ejection feedback (dot- 
ted), ejection feedback + dust (dashed). Lower panel: ACDM, retention feedback (solid), retention 
feedback + dust (dashed). 



in clusters relative to the field and as a result, ^(r) is steep on small scales. 
6.5 Pairwise Velocity Moments 

Figure 13 shows the first (^12) and the second moment (0-12) of the one-dimensional radial 
pairwise velocity distribution of galaxies and dark matter particles as a function of pair sep- 
aration for our best-fit models. Bold lines are the moments for the dark matter particles. 
Solid and dotted fines are the moments for galaxies brighter than Mb — —17.5-1- 51og/i and 
Mb — —18.5-1- 51og/i, respectively. Brighter cutoffs yield similar results. Galaxies do not 
show any significant velocity bias on any non-linear scale regardless of their luminosity. This 
result is not trivial because galaxies are not a random subsample of dark matter particles. 
This result is consistent with the weak density bias and the lack of luminosity segregation 
described in Sect. 6.6. 

The pairwise velocity dispersion profile ai2{r) is sensitive to the presence of clusters within 
the sample and therefore to the total volume of the sample. Our simulation box has a volume 
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of the order of lO^h ^ Mpc^ which is comparable to the volume ~ 7 x W^h ^ Mpc^ within 
cz = 12000 km s^^ of the northern slice of the Center for Astrophysics Redshift Surveys 
(hereafter CfA2N, de Lapparent et al 1991; Geller & Huchra 1989; Huchra et al 1990; Huchra 
et al 1995). Filled squares in figure 13 show ai2{r) for the CfA2N shoe. Both rCDM and 
ACDM models are in remarkably good agreement with the CfA2N slice. 

Lacking the full three-dimensional information, we compute the CfA2N values by fitting 
the redshift space correlation function ^{rp,7r) with the convolution model (Fisher 1995): we 
model C(rp, tt) as the convolution of the real space correlation function with an exponential 
pairwise velocity distribution (Davis & Peebles 1983; details of this procedure are given in 
Marzke et al 1995). In order to apply this procedure, we need to assume a model for the mean 
streaming velocity Vi2{r). We adopt the similarity solution suggested by Davis & Peebles 
(1983) Vi2{r) = r/[l + (r/ro)^] where tq = 5.83h~^ Mpc is the CfA2N galaxy-galaxy correlation 
length. Assuming Vi2{r) = usually yields lower values of ai2 (see e.g. Marzke et al 1995; 
Somerville, Primack & Nolthenius 1997). Note that our simulations show Vi2 ~ 200 — 300 km 
s~^ at relative separations r ~ lh~^ Mpc and do not support the assumption of a null Vi2{r). 

In future papers, we use mock redshift surveys to show that this procedure yields pairwise 
velocity dispersion profiles 0"i2(r) in reasonably good agreement with the real profile. It is 
therefore encouraging that the a^ir) profiles of our best- fit models match the CfA2N result. 

6.6 Galaxy Bias 

An important issue we wish to address is that of galaxy "bias", i.e. whether the clustering 
amplitude of the galaxies in our models differs from that of the underlying dark matter. In 
figure 14, we compare ^{r) for galaxies of different luminosities, morphological types and 
colours with ^(r) for the dark matter. For rCDM, we show the model with ejection feedback 
and dust extinction and for ACDM, we show the model with retention feedback and dust 
extinction. These are the "best fit" models in both cases. . 

On large scales, galaxies of all luminosities are unbiased in both models. On small scales, 
the rCDM galaxies are shghtly more clustered than the dark matter, whereas the ACDM 
galaxies are a factor ~ 2 — 3 less clustered. Neither model displays any luminosity segregation; 
galaxies of all magnitudes cluster in the same way. As discussed by Kauffmann, Nusser & 
Steinmetz (1997), positive bias is obtained for galaxies in halos more massive than M^{z = 0), 
the mass of the typical collapsed object at the present day. For cluster-normalized CDM 
models, is large (~ lO^^M©). Our bright normalization also causes luminous galaxies to 
be placed in relatively low mass halos. This is why we do not see the luminosity-dependent 
bias discussed by Kauffmann, Nusser & Steinmetz (1997). 

On the other hand, red galaxies and to a lesser extent early-type galaxies are more clustered 
than the underlying population, particularly on small scales. These galaxies occur preferen- 
tially in the cores of groups and clusters (see figures 4 and 5). Conversely, star- forming galaxies 
are less clustered on small scales because they are predominantly central halo galaxies and 
are thus spatially exclusive. A similar result can be seen in the earlier work of Cen & Ostriker 
(1992). 

The relative bias between different types of galaxies in the Southern Sky Redshift Survey 
has recently been analyzed in a paper by Willmer, Da Costa & Pellegrini (1998). These authors 
find that galaxies with magnitudes below ~ do not exhibit any luminosity-dependent 
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Figure 13: Mean streaming velocity Vi2 (lower curves) and pairwise velocity dispersion 
(Ti2 (upper curves) at different relative separations for our best-fit models. Bold, solid, 
and dotted lines are the moments for the dark matter particles, galaxies brighter than 
Mb — —17.5 + 51og/i and Mb — —18.5 + 51og/i, respectively. Sohd squares are the pairwise 
velocity dispersions for the CfA2N redshift survey. Error bars are computed with a bootstrap 
resampling method. 
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Figure 14: Galaxy bias as a function of luminosity, type, colour and star formation rate. The left 
panels axe for the rCDM simulation, the right panels are for the ACDM simulation. Upper panels: 
The thick solid line is the dark matter correlation function. Thin solid lines show ^(r) for different 
B-band luminosities: 5 < -19 (solid), B < -20 (dotted), B < -20.5 (short-dashed), B < -21 
(long-dashed). Lower panels: Thick solid line is the dark matter correlation function. Thin solid 
line corresponds to early-type galaxies, dashed line to red {B — V > 0.8) galaxies, and dotted line 
to star forming {SFR > 2Mq yv~^) galaxies. 

bias, but at magnitudes brighter than L*, galaxies are significantly more strongly clustered 
(6(L > L^)/h{L^) = 1.5). As discussed above, our models do not exhibit this luminosity- 
dependent bias, perhaps indicating that our adopted normalization is too bright. Willmer, 
Da Costa & Pellegrini also find that both early-type galaxies and red galaxies cluster more 
strongly than the galaxy population as a whole. The bias is stronger on small scales (r < 4 
Mpc) than on larger scales. In addition, galaxies with colours characteristic of old stellar 
populations exhibit more bias than galaxies with early-type morphologies. These results are 
in good qualitative agreement with the results shown in figure 14 and the earlier results of 
Kauffmann, Nusser & Steinmetz (1997). 
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6.7 Mass-to-light Ratios of Clusters 

Clusters of galaxies are the largest collapsed objects in the universe and have often been used 
as a means of estimating the cosmological density parameter Q. The standard procedure is as 
follows. The mass of a cluster within some radius is measured using techniques such as galaxy 
kinematics, X-ray profiles or gravitational lensing. This is compared to the luminosity within 
the same radius to calculate a mass-to-light ratio. The mean mass density of the universe is 
then estimated by multiplying its mean luminosity density by this number. The validity of 
this technique depends on two questionable assumptions: 

1. Luminosity is "conserved" when field galaxies fall into a cluster. In practice, we know 
that cluster galaxies are predominantly ellipticals with no ongoing star formation, whereas 
field galaxies are predominantly spirals in which stars are forming at several solar masses 
per year. The stellar mass-to-light ratios of field and cluster galaxies thus differ strongly 
at optical and UV wavelengths. 

2. The efficiency of galaxy formation is the same in all environments. 

Most studies of cluster mass-to-hght ratios give values corresponding to ~ 0.2 (see for 
example Carlberg et al 1996). If a high-density cosmology is to be consistent with these 
observations, we require the mass-to-light ratios of clusters to be significantly lower than the 
mass-to-light ratio of the Universe as a whole. This is often also referred to as "bias" in the 
literature, though this kind of bias is quite different to what we measure when we compare 
the correlation functions of galaxies and dark matter. 

Figure 15 shows the mass-to-light ratios of our simulated halos divided by the mass-to-light 
ratio of the simulation as a whole. The results are shown as a function of the virial mass of 
the halo and for both B and I-bands. It is interesting that in the Q = 1 rCDM model, the 
mass-to-hght ratios of clusters are significantly lower than the global value. The factor (~ 2) 
that we obtain is somewhat too small, however, to bring the model into good agreement with 
the observations. Rich clusters in the ACDM model have mass-to-light ratios only 10-20 % 
smaller than the global value. It is also interesting that groups of galaxies are predicted to 
have lower mass-to-light ratios than clusters in both models. This agrees with the observed 
trends. 

6.8 Colour Distributions 

We now compare the B — V colour distributions of our model galaxies with colours of bright 
nearby galaxies in the RC3 catalogue (De Vaucouleurs et al 1991). We have selected 1580 
galaxies with Mb < —19.0 — 51og/i, which have both morphological classifications and mag- 
nitudes in the Johnson B and l^-bands. Of these galaxies, 653 are classified as early-type. 
Note that the fraction of early-type galaxies is rather high, which may reflect the fact that the 
sample is dominated by galaxies in the Virgo cluster. In our simulation ~ 25% of galaxies to 
this magnitude limit have M{B)i,uige — M{B) total < 1 ^lag and so are classified as early types. 

In figure 16, we compare the colour distributions of galaxies in the rCDM model with that 
of the RC3 sample. The left panel shows the model with ejection feedback, but without any 
dust extinction. The right panel shows the same model with dust extinction. Without dust 
extinction, the model colour distribution peaks at B — V ^ 0.6 in contrast to the observed 
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Figure 15: The mass-to- light ratios of halos divided by the mass-to-hght ratio of the Universe 
as a whole. The results are shown as a function of the virial mass of the halo. 
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distribution, which is relatively flat from B — V ^ 0.4 to B — V ^ 1. In particular, early- type 
galaxies appear to be significantly too blue in this model. With the inclusion of dust, the 
color distribution broadens significantly. Many of the field ellipticals on our models have a 
lower-luminosity, star-forming disk component. The inclusion of dust in the models reddens 
these objects and brings the colour distribution of early-type galaxies into better agreement 
with the observations. 

The colour distributions of galaxies in the ACDM models are shown in figure 17. Because 
galaxies form earlier in this cosmology, the colours tend to be redder on average, but the 
effect is small. Again, the effect of dust extinction is to broaden the colour distributions and 
to bring them into better agreement with the observations. 

6.9 Star Formation Rate Functions 

In addition to the luminosity function and the colour distributions of galaxies, observational 
estimates of the number density of galaxies as a function of star formation rate provide 
an important constraint on galaxy formation theories. Gallego et al (1995) have analyzed 
a complete sample of emission-line galaxies and have computed an Ha luminosity function. 
Assuming a Scalo (1986) IMF and using Case B recombination theory to predict the luminosity 
of the Ha emission line for a given rate of star formation, they derive an integrated star 
formation rate density of the local Universe of 0.013 ± O.OOlMgyr"^ Mpc ~^ (for Hq — 50 km 
s"^ Mpc"^). 

In figure 18, we compare the Ha luminosity functions of our models with the results of 
Gallego et al. We have used the same transformation as these authors to convert from star 
formation rate to Ha luminosity, 

L{Ha) ^ 9.40 x 10^° ^ ^^-i ergs s^^ (17) 

The results of the rCDM model agree better with the observations than those of the ACDM 
model. For rCDM, we obtain an integrated SFR density of 0.017 M© yr~^ Mpc~^. For ACDM 
the present-day SFR density is a factor ~ 4 too low. From figure 18, we see that there there 
is a lack of galaxies forming stars at rates between 1 and 10 Mq yr~^. These are the same 
galaxies that are missing at the "knee" of the B-band luminsity function in figure 9. 

7 Discussion and Conclusions 

None of the models analyzed in this paper provides a completely satisfactory fit to all the 
observed properties of galaxies at 2; = 0. If supernova feedback is weak, the = 1 rCDM 
model produces too high a total luminosity density. Most of the excess luminosity is in the 
form of star-forming "field" galaxies. As a result, the galaxy correlation function turns over 
on scales below a few Mpc. The observed two-point correlation function, on the other hand, 
is well represented by a single power law of slope —1.7 down to scales approaching 10 kpc 
(Baugh 1996). If one assumes instead that supernova feedback is able to eject gas out of a 
halo, so that it becomes unavailable for further cooling and star formation until much later 
times, the number of star-forming galaxies in low-mass field halos is significantly lowered and 
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Figure 16: The B — V colour distribution of galaxies brighter than B = —20.5 in the rCDM 
simulation. The left panel shows a model without dust and the right panel a model with dust. 
Dotted lines show the model distributions. Solid hues show the distribution oi B < —20.5 
galaxies with both colours and morphologies in the RC3 catalogue. 
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Figure 17: The B — V colour distribution of galaxies brighter than B — —19.7 in the ACDM 
simulation. The left panel shows a model without dust and the right panel a model with dust. 
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Figure 18: The Ha luminosity function of galaxies in the simulation compared to the 
Schechter fit derived by Gallego et al (1995). 

this improves the fit to the observed luminosity and correlation functions. Even so, the model 
still produces an excess of very bright galaxies in groups and clusters. The Q = 0.3 ACDM 
model has the opposite problem: it produces too few galaxies, even if feedback is weak and 
no gas leaves the halo. The two point correlation function is also too steep on scales below a 
few Mpc. 

One obvious inference we can draw from these results is that a CDM model with density 
parameter somewhere between 0.3 and 1 may well do better in matching the observations. 
Indeed, in recent work, Somerville & Primack (1998) find that a model with Q = 0.5 provides 
the best fit. 

Care must be exercised, however, before drawing strong conclusions about cosmological 
parameters. Many of our difficulties in fitting the observed shape of the galaxy luminosity 
function come about because we have chosen a very bright normalization for the models. 
This is because we force the models to fit the zero-point of the Giovanelli et al I-band TuUy- 
Fisher relation at a circular velocity of Vc — 220 km s~^. We have also assumed that the 
circular velocities of central disk galaxies are the same as those of their surrounding halos. In 
practice, the relation between the circular velocities of disk and halo is not straightforward 
and requires a detailed model for the internal structure of the halo and of the galaxy within it. 
More detailed calculations of disk formation in CDM halos predict that the circular velocity of 
the disk is generally higher than that of the surrounding halo (Mo, Mao & White 1998). This 
results in a higher total luminosity density for the same TuUy-Fisher normalization, and would 
improve the fit of the ACDM model while making the rCDM model even more discrepant. 
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One might also question whether normahzing the models to match the Tully-Fishcr relation 
is sensible, since the physical origin of this relation is still uncertain and depends on structural 
properties of the galaxies which we are not attempting to model. We have argued that the 
distribution of K-band luminosities should provide a robust test of galaxy formation theory, 
since the K-band luminosity of a galaxy depends primarily on its stellar mass and not on its 
instantaneous star formation rate or its dust content. With the advent of wide-field K-band 
galaxy surveys such as the Two Micron All Sky Survey (2MASS, Kleinmann et al 1994), 
normalizing the models to match the observed K-band luminosity function may be the best 
choice for the future. 

Another effect we have investigated is that of dust extinction in star-forming galaxies. If 
dust is ignored, the Q = 1 rCDM model produces too many blue galaxies compared to the 
observations. Galaxy formation occurs rather late in this model, so the stellar populations of 
galaxies are younger than in the ACDM model. However, the inclusion of a simple, empirically- 
motivated dust extinction scheme makes this problem go away, and the colour distributions 
of galaxies in the rCDM and ACDM models are then virtually indistinguishable. Since dust 
extinction only affects star-forming field galaxies, models with dust produce somewhat steeper 
and higher amplitude correlation functions. Finally, dust lowers the predicted number of 
galaxies at the bright end of the B-band luminosity function. Our models tend to produce a 
power-law "tail" of high-luminosity galaxies, rather than a sharp exponential cutoff. Although 
dust can reduce this tail in the B-band, the problem persists in the K-band. 

Finally, although we have tried to illustrate the effect of different schemes for star for- 
mation, feedback and dust, we cannot be sure that we have taken into account all the key 
physical processes that may seriously affect our results. For example, we have proposed that 
very efficient feedback can stop galaxies forming in low-mass halos, but alternative schemes 
may be possible. For example, Jimenez et al (1997) have proposed that the gaseous disks 
that form in low Vc halos are not subject to the instabilities that cause star formation and 
are thus not seen in the optical or UV, although their emission at 21 cm wavelengths should 
be detectable. The only way to make progress is through observational studies of how star 
formation rates, cold gas and dust contents, and the shells and bubbles seen around star- 
forming regions, vary with galactic mass, metallicity, morphology and environment. This is 
a massive and long-term undertaking, but it seems necessary if we are to build a picture of 
galaxy formation that is more than just schematic. 

We conclude that we cannot reliably constrain the values of cosmological parameters using 
the properties of the galaxy distribution at z = because our results are strongly infiuenced 
by our adopted normalization and the way in which we choose to parametrize star formation 
and supernova feedback. The main effects of changing the cosmological parameters in our set 
of cluster-normalized CDM models are 

1. to shift the epoch when halos of a given mass form 

2. to change the number density of galactic- mass halos at the present day. 

The adopted normalization and the recipes for feedback and star formation affect the efficiency 
with which stars form as a function of redshift and galaxy mass. It is thus not surprising that 
there is considerable freedom to "tune" our results to fit the observations. Nevertheless it 
is interesting that very different star formation and feedback schemes are required to come 
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close to reproducing the observations in the two cosmologies we have explored. For = 1, 
we require that feedback be very efficient at removing gas both from galaxies and from their 
surrounding halos. Many lower-mass halos are "dark" in this model, containing a faint, 
possibly low-surface brightness object. In a low-density Universe, feedback must be inefficient 
and every halo must contain a star-forming galaxy if we are to match the observations. This 
gives us confidence that if the cosmological parameters of the Universe can be determined 
using techniques such as the analysis of small-scale fiuctuations in the microwave background 
(e.g. Bond, Efstathiou & Tegmark 1997) or the light-curves of high-redshift supernovae (e.g. 
Perlmutter et al 1997), the methodology introduced in this paper will become valuable for 
understanding and constraining the detailed physical processes operating within galaxies. 

It is also encouraging that our best-fit models come reasonably close to fitting many 
different observed measures of the galaxy distribution at the same tim,e. The models that 
give the best fits to the B-band luminosity function, also give the best fits to the K-band 
luminosity function, the correlation function, the color distribution and the star formation 
rate function. In Paper II we demonstrate that these same models provide a reasonably 
good match to the observed properties of galaxy groups. This gives us confidence that the 
fundamental theoretical framework, which determines how galaxy properties scale with halo 
mass, may indeed be correct. 
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Table 1: Cosmological parameters of the GIF models. Qq and JIa are the density parameters 
for matter and for the cosmological constant, h is the Hubble parameter, erg is the rms of the 
density field fluctuations in spheres of radius 8 Mpc, and F is the shape parameter of the 
power spectrum. Also given is the size of the cosmological simulation box. 



Model 






h 


0-8 


F 


Box Size 














[Mpc/h] 


SCDM 


1.0 


0.0 


0.5 


0.60 


0.50 


85 


rCDM 


1.0 


0.0 


0.5 


0.60 


0.21 


85 


ACDM 


0.3 


0.7 


0.7 


0.90 


0.21 


141 


OCDM 


0.3 


0.0 


0.7 


0.85 


0.21 


141 
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